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Abstract The four-point integral of the minimal super Liouville gravity on the sphere 
is evaluated numerically. The integration procedure is based on the effective elliptic param- 
eterization of the moduli space. The analysis is performed for a few different gravitational 
four-point amplitudes. The results agree with the analytic results recently obtained using 
the Higher super Liouville equations of motion. 

1 Introduction 

The continuum formulation of the noncritical string theory is equivalent to 2D quantum grav- 
ity coupled to some critical matter, i.e., the matter described by a conformal field theory 
Aic- Simple reaction of conformal theories to the scaling of the metric leads to the universal 
form of the effective action of the generated gravity, which is called the Liouville gravity 
(LG) [1]. Because of the peculiarities of two-dimensional metric geometry, many technical 
simplifications immediately come into play. Thus, LG is perhaps the simplest example of 
quantum gravity, but it nevertheless shares the same basic questions of interpretation and 
can hence be considered useful and worth studying. The problem of choosing observables 
correctly and the problem of calculating the corresponding correlation functions are of pri- 
mary importance in any quantum theory. The field of LG has experienced considerable 
progress in recent years. Recently discovered higher equations of motion (HEM) [2] in LET 
have allowed reaching the four-point level in calculating the correlation functions in LG [3, 4]. 
The results were tested against the calculation in the framework of the relatively indepen- 
dent approach to 2D quantum gravity usually called the matrix models (see, e.g., [5] and the 
references therein). Moreover, a deeper understanding of the correspondence between these 
two approaches was achieved based on these results [6] , although the complete picture of the 
relations between these techniques is still missing. 

In the context of string theory applications, the construction of 2D quantum gravity in 
superspace is one of the most interesting questions. The first possible generalization is A^=l 
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supersymmetry. Here, essential progress was also achieved recently. The study was motivated 
by a series of works concerned with super LFT in which the profound understanding of 
the properties of the conformal blocks (which are the basic elements of any CFT [7]) was 
achieved. This allowed "taking off" the necessity of treating the correlations involving only 
special degenerate [8] excitations. Another important result was the discovery of the HEM 
in super LFT [9] . All these results have served as a good starting point for a more profound 
study of super Liouvillc gravity (SLG). In [10], the structure of the physical fields in the 
Neveu-Schwarz (NS) sector of SLG was clarified, and the general expression for the n- 
point correlation numbers on the sphere in terms of integrals over moduli space was written 
exphcitly. Then the super HEM in the A^=l super LFT and the analysis of the structure of 
the super ground ring (and its logarithmic counterpart) allowed deriving the explicit analytic 
expressions for the four-point correlators. Because little is now known for the supersymmetric 
matrix model, there are no independent results analogous to those obtained in [10]. In this 
situation, more checks of the validity of these results would seem desirable. This paper is 
devoted to directly calculating the four-point correlation numbers in SMLG. 

The paper is organized as follows. To make the presentation self-contained, we collect all 
necessary information related to the subject in Sec. 2 and Sec. 3. The remaining part of the 
paper deals with evaluating the four-point integrals numerically. In Sec. 4, we consider two 
examples of the four-point integral. We reduce the expressions to a number of integrals over 
the fundamental region of the modular group, and the integrands represent the products of 
the various correlation functions both in the matter and in the Liouville sectors. The elliptic 
transformation is the subject of Sec. 5. The numerical results arc presented in Sec. 6. Some 
useful details omitted from the main text are given in Appendix A. 

2 Minimal Super Liouville Gravity 

In the framework of the so-called DDK approach [11-13], SLG is represented as a tensor 
product of superconformal matter (SM), super Liouville, and super ghost systems 



^SLG — ^SM + ^SL + ^SG 



(1) 



with the interaction via the relation for the central charge parameters 



CSM + CSL + CSG = 



(2) 



and also due to the construction of the physical fields. 
The superconformal algebra is 

[L„, Lm] ^ {n- m)Ln+m n)8n, 




(3) 
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where ^ 

r, s e Z + - for the NS sector, 

r,s for the R sector. 

The SLFT central charge is 

csL = l + 2g^ (4) 

where the "background charge" parameter Q is related to the SLFT basic quantum param- 
eter h 

Q^b-^ + b (5) 

The fields belong to the highest-weight representations of the superconformal algebra. The 
basic fields of interest in this paper belong to the primary supermultiplet {Va, Ya, Ya, Wg) 
with the bottom component T4 having the conformal dimension 

= (6) 

where a is a continuous (complex) parameter, and the other components of the primary 
supermultiplet being 

Ya — G^i/2^a, Ya = (5^1/2^, Wa = G^l/2G^l/2^a- (7) 

Here and hereafter, we use the superscripts M, L, and G to specify the matter, Liouville, and 
ghost sectors of the superconformal generators. The generators without sector superscripts 
are related to the total super Virasoro algebra. 

At certain special values of the parameter a — am,m one singular vector appears at the 
level mn/2 in the Verma module over Kr„,n — Vm,n [8]- Here, 

where (m, n) is a pair of positive integers {m — n & 2Z) and 

A„,„ = =^. (9) 

The basic super Liouville operator product expansion [14] (for the sake of brevity we 
write A = ^Q/2+iP and A, = A„.) 

Va,{x)VaM= (10) 
/' ^ (XX)— (Cg/^:- [V,/2..A0)L + «r [%2..P(0)]„J 

This OPE is continuous and involves integration over the "momentum" P. In (10), [V^lg^oo 
denotes the contribution of the primary field Vp and its "even" and "odd" superconformal 
descendants to the operator product expansion. As usual, the prime on the integral indicates 
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possible discrete terms; in this study, we consider only the region b where such extra terms 

do not appear and the integral can be understood literally. All other OPEs of two arbitrary 
local fields in SLFT can be derived from (10). The basic structure constants and 
in (10) were evaluated using the bootstrap technique in [15, 16] and have the exphcit 
form (here a denotes Oi + 02 + 03) 



13 



as 
aia2 



{Q-a)/b 



TR(6)TNs(2ai)TNs(2a2)TNs(2a; 



2TNs(a. — (5)TNs(ai+2-3)'^Ns(a.2+3-l)'^Ns(fl3+l-2) ' 



{Q-a)/b 



iTR(6)TNs(2ai)TNs(2a2)TNs(2a3) 



(11) 



TR(a - (9)TR(ai+2-3)TR(a2+3-i)TR(a3+i_2) ' 
where we use the convenient notation in [17] for the special functions 



X + b 



X + b 



(12) 



expressed in terms of the "upsilon" function T^,, which is the standard element in the Liouville 
field theory (see [18,19]). 

Because of central charge balance condition (2), the central charge of the matter sector 
is given in terms of the same basic parameter b: 



csM = 1 - 2gf^ 



(13) 



where q — b ^ — b. We let ($q,, Xa,Xa, ^a) denote the primary multiplet in the matter sector 
with the dimension of the bottom component being given by 



,M 



a{q — a) 



(14) 



The super ghost system (see, e.g., [20-22]) is described by the free super conformal field 
theory with the central charge csg = — 10. The fermionic part of the SG system involves two 

anticommuting fields {b,c) of spins (2, —1), and the bosonic part involves two bosonic fields 
(/5,7) of spins (3/2, —1/2). The formal fields (see [10]) of the form ^(7(0)) of dimension 1/2 
are essential in constructing the gravitational amplitudes. 



3 Physical Fields and the Correlation Numbers 

The physical fields form a space of cohomology classes with respect to the nilpotent BRST 
charge Q, 



-M+L 



7- 



(15) 



4 



In this work, we deal with the correlators of physical fields of the two types 

Waiz, z) = Ua{z, z) ■ c(^)c(^) • 5{^{z))5{^{z)) , (16) 

and 

W„(^, -z) = [c'^tk + ^^-1/2) (g^J^ + ^G^_v2)ua(^, -z) ■ c{z)c{z), (17) 

where 

\]a{z,z)^^a-b{z,z)V,{z,z) (18) 

Here the parameter a can take generic values. The general form of the n-point correlation 
numbers on the sphere for these observables [10] is 

4(ai,--- ,an) =ny" ^i'^i(|G_v2G'_v2U„,(^,)W„3(z3)W„,(^2)W„,(zi)y (19) 

An additional "discrete" physical state arises when the representation in the matter sector 
is degenerate, 

Om,„(^, Z) = Hm,nHm,n^m,n{z, z)Vm,n{z, z). (20) 

The operators Hm,n are composed of the super Virasoro generators and arc defined uniquely 
modulo exact terms. Moreover, if we introduce the logarithmic counterparts of the discrete 
states Om,n, 

= Hm,nHrn,n^m,nV:,^n, (21) 

then we have the important relations [10] 

QQO:„_„ = B^^^Wm,-n (22) 

and 

G_v2G-i/2U„,_„ = B-]JdO'^^^ mod Q, (23) 

where „ are the coefficients arising in the higher equations of motion of SLFT [9]. For 
four points, relation (23) allows reducing the moduli integral in general expression (19) to 
boundary integrals if one of the fields is degenerate, i.e., = am,-n- The explicit result is 
(see [10]) 

h{am,-n:ai,a2,a3) = KN{am,-n)lY^ qi'^'''\ai) + 2mnXm,n\YlN{ai), (24) 

^ j=l r,se(m,n) ^ i=l 



where 



g(™'")(a) = \a- Xr,s - g/2|Re - Xm,n (25) 
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and the fusion set is (m, n) — {1 — m : 2 : m — 1, 1 — n : 2 : n — 1} . The coefficient is 



and the "leg" factors are 

N{a)= [7r/x7(^2+ 2" 

4 Direct Calculation 







2+b-2 


"r^ 1" 
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(26) 



-a/b 



1/2 



(27) 



Here, we verify analytic result (24). The space of parameters (ai, a2, 03, and also b) is rather 
big to present a comprehensive analysis. In what follows, we focus on the two examples 
where one of the fields is either or W2f,. Moreover, we restrict ourself to considering 
the most symmetric situation of four identical fields ^^{a) — Ii{a, a, a, a). In the four-point 
case (19) reduces to 



(28) 



X4(a)= j rf'z(G'_i/2G-i/2U„(;2)W„(0)W„(l)W„(oo)). 
The analytic results following from expression (24) are 



where 



For the second integral, we have 



I,{b) = lN\b)Y.^'^'\b), 



E(^'1)(6) = |26' - 1|. 



and 



X^{2b) = '^N\2b)i:^''''\b) 



\5b^ - 11 + 136^ - 11 + 16^ -11-36^-1 



(29) 



(30) 



(31) 



(32) 



Let us consider the integral I^ib). Taking into account that we deal with the unit oper- 
ators in the matter sector in this case, we have 



G_i/2G_i/2Ufe = Wh-, 
W, = Hcc<5(7)5(7), 



(33) 
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Taking the explicit form of the correlation functions in the ghost sector into account, 



(C(0)C(1)) = 0, 
(C(0)C(l)C(oo)) = 1, (34) 
(5(7(0))5(7(1))) = 1, 

we conclude that the only nonzero contribution comes from the term in W^, which is pro- 
portional to cc, 

2-4(6) = j d^z{Wh{z)V,{Q)W,{l)V,{oo)). (35) 

In the same way, we have 

W2b^^bV2bCc5{^)5{^), (36) 

for the second integral, and taking (34) into account, we obtain 

U2b)^ J rf2;,|^(v]/,(^)$,(0)*,(l)$,(oo))(y26(^)^26(0)y26(l)^26(oo)), 

+ (X6(;^)$6(0)X6(l)$6(oo))(F26(^)V^25(0)F25(l)V^2fe(oo)), (37) 
+ (Xb(^)$6(0)X6(l)$b(oo))(F2b(^)V^2b(0)F2b(l)V^26(oo)), 

+ ($;,(Z)$5(0)$6(1)$;,(00))(W^26(-2)^26(0)W^2(,(1)^2(,(00)), 

We now use the symmetry of the integrals under modular transformations to reduce 
the integration from the whole complex plane to the fundamental domain. The modular 
subgroup of projective transformations divides the complex plane into six regions. The 
fundamental region is defined as G ={Rex < 1/2; |1 — a:| < 1}. The other five regions 
are mapped to the fundamental one using one of the transformations A, B, AB, BA, ABA, 
where A: z ^ 1/z and B: z 1 — z. Combining the projective transformations of the fields 
and the corresponding change of variables in the integrals, we reduce the integration to the 
fundamental region. We note that the Jacobian of the transformation exactly cancels the 
transformation of the fields because their total conformal dimension is 1. Then, 



2-4(6) = 2 / d'z{ {Wb(z)Vb(0)Wb(l)Vb(oo)) + {Wb(z)Vb(0)Vb(l)Wb(oo))+ 

+ {Wi,{z)Wh{0)Vb{l)Vb{oo))] , (38) 
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where the factor 2 in front of the integral takes the equivalent projective images into account. 
The expression for the second integral is rather bulky: 



Zi(26) ^2 d^z 
Jg 



{^,{z)MQ)^b{l)Moo)){V2b{z)V2b{0)V2b{l)V2b{oO 

+{M^)^b(0)M'^)Moo)){V2b{z)V2b{0)V2b(l)V2b(oo 

+ {^b{z)M0)M^)^b{oo)){V2h{z)V2b{0)V2b{l)V2b{oo)) 



+ \^{Xbiz)MQ)Xbil)M^)){y2biz)V2i>iO)Y2i>il)V2bioo 

+ (X6WX6(0)*6(l)*6(oo))(F26(^)F2ft(0)^6(l)V^26(oO 
+ (X6(^)$6(0)$6(l)x6(oo))(y26(-2)^26(0)y26(l)?'26(oo)) 



+ \^{Xb(z)Mmb{i)M00)) {Y2b(z)V2b{0)Y2b{l)V2b(oO 

+ {Xb{z)Xb{Q)M^)Moo)){Y2b{z)Y2i>{0)V2b{l)V2b{cK> 

+{Uz)Mo)Mmb{<^)){y2b{z)V2b{o)V2b{i)Y2b{<x,)) 



+ \^{^b{z)MQ)M^)Moo)){W2b{z)V2b{0)W2b{l)V2b{oO 
+ {^b(z)M^)M'^)Moo)){W2b{z)W2b(0)V2b(l)V2b(oo) 

+ {M^)^b{0)M^)Moo)){W2b{z)V2b{0)V2b{l)W2b{oo)) 



(39) 



We now use the conformal block decomposition of the correlation functions. It is useful 
to introduce a compact notation. For a while, we omit some arguments that are easily 
reconstructed in the final expressions. In the matter sector, 



($(.)$(0)$(l)$(oo)) = c,|Af(^)r, 

{z)m^ mioo)) = c,\Ai\z)\\ 
{z)^ mmioc)) = c,\A^;\z)\^ 
{nz)mmnoo))^c,\A'-^\z)\'. 



(40) 



Here, the index k — +, 0, — corresponds to the three channels in the degenerate OPE of the 
field $6 (and also of its super partners), and we assume summation with respect to k. The 
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coefficients Ck are related to the basic structure constants (see [10]): 



. r^(h^ 7(1/2 + &V2)7(-l/2 + 5&V2) 

+ 7(-l/2 + 362/2)7(1/2 + 3672)' ^ ^ 

Co = -C',{b) = -1, (42) 
^ - (,) - 7(2b^)7(-l/2 + 672)7^(1/2 + &V2) 

-U b^^^{b^)j{-l + b^)j{-l/2 + 3by2)' ^ ' 

In (40), A^"^ denotes the conformal blocks appearing in the k channel for the given corre- 
lation function. Here and hereafter, the normalization is chosen such that all but the basic 
combinations are absorbed inside the conformal blocks. In Appendix A, we recapitulate 
some details and explicit constructions concerning conformal blocks. In the Liouville sector, 

{V{z)vm{i)v{<^))^n j ^n{p)\Bl'\p,z)\', 
{W{z)v{o)w{i)v{oo))^n I ^n{P)\Bl'\p,z)f, 

/fi p 
—ri{P)\BP{P,z)\', (44) 

{W{z)V{0)V{l)W{oo))^n J '^r,{P)\Bl'\p,z)\'. 

Because the correlation functions do not contain the degenerate fields in the Liouville sector, 
the index / here assumes summation of the "even" and of the "odd" conformal blocks in 
accordance with general OPE (10). Again, the normalization leaves only two basic combi- 
nations outside the conformal blocks, 

^n(P) = Cgr^^Cgr^^, (45) 

where we separate the factor TZ for convenience (it is independent of P) and the parameter 
a is related to the external conformal dimension (i. e. either b or 26). All other correlation 
functions in (39) are already not independent. They are expressed in terms of the same 
ingredients as in (40) and (44). For example, 

(x(^)$(0)x(l)$(oo)) = c,A«(z)4°)(z), 

(X(^)x(0)$(l)$(oo)) = c,A^^\z)Af\z), (46) 
{x{z)mmx{oc)) = c,A^^\z)A^^\z), 
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and 



{Y{z)V{0)Y{l)V{oo))^n J —ri{P)Bl^'{P,z)Bl''{P,z), 
{Y{z)Y{0)V{l)VM)=n I ^ri{P)Bl'\P,z)Bl'\P,z), (47) 



dP 



{Y{z)V{0)V{l)Y{oc))^n J "—n{P)Bl'\P,z)Bl^>{P,z). 

The remaining six correlation functions not written explicitly are obtained from (46) and (47) 
by complex conjugation. Using the introduced notation, we can rewrite the integrals under 
consideration in the compact forms 



j(3), 



Bi'\P,z)\' + \Bl-^>{P,z){' + \Br>{P, 



?(2) 



(3), 



(48) 



and 



U2b)^2n I £z I ^j;c,rKP)[|A«(^)sJ°)(P,^)+Af(^)S«(P,^)r 
"^^ k,i L 

+ \A^^\z)b\'\p, z) + Af{z)Bf\P, z)\' + \Af{z)Bf\P, z) + Af{z)Bf\P, z){^ 



(49) 



Bulky expression (39) has a remarkably compact and clear structure in terms of the conformal 
blocks. We again note that in (48) and (49), we respectively assume the different values of 
external conformal dimensions in the Liouville sector and A26. 

5 The Modular Integral 

It turns out efficient [23] to use eUiptic transformations in the integration. We use the 
standard map 



T — l- 



K{l-z) 



K{z) ' 

where the complete elliptic integral of the first kind is 



(50) 



Kiz) 



1 T— 



(51) 



and y"^ = t{l — t){l — zt) . It can be verified that 



dz — TTz{l — z)9^{q)dT, 



(52) 
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where 



and 



n=— oo 

Integral (48) becomes 

^(l-^)^3^(g)sf)(P,^)|^dV+ / \z{l - z)9t{q)Bl'\P, z)\'d' 

J F 



where F={|t|>1; |ReT|<l/2}. Similarly, for (49), we have 



no 47r 



/■ |^(l-;^)^3^(g)(4^)(;.)i?;°)(P,;^)+4°)(^)s(^n^,^))r^V 
+ /■ |z(l-z)^3^(g)(4^)(^)i?;°H^,^) + 4°^(^)i?fH^,^))r^'r 

«/ F 

/■ 1^(1 - z)el{q){Af\z)Bi'\P,z)+Af\z)Bi'\P,z))\'d'T 

J F 



+ 



We now define the conformal blocks more explicitly, 



J^^{z) 


= i^o^(0,^), 




= F^,{b, 


z), 


Afiz) 


= F^{2b, 


z), 


Ay}{z) 


= i^n(0,^), 




= F^oib, 


z), 


Ai\z) 


= F^,{2b, 


z), 


J^^{z) 






= F^Ab, 


z), 


A?\z) 


= F^{2b, 


z), 


J^^{z) 


= F^iiO,z), 




= F^oib, 


z), 


Af{z) 


= Fr^i2b, 


z). 



(53) 



(54) 



(55) 



(56) 



(57) 



Here, the first argument of the symmetric conformal blocks defines the internal conformal 
dimension. The first lower index corresponds to one of the four basic types of conformal 
blocks wc consider with respect to the set of external fields (sec Appendix A); the second 
index is if the corresponding block with the given internal conformal dimension is "even" 
and 1 if it is "odd" . In the Liouville sector, 



bT (^) = ^o'o {P, z) , < iz) = i^o'i {P, z) , 

Bi,'\z) = F^i(P, z), Bi'\z) = Fi,(P, z), 

B'i\z) = F^,{P, z\ Bf{z) = Fjl(P, z\ 

Bf\z) = Fl,{P, z\ Bf{z) = P3VP, z). 



(58) 
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For the first integral, the comphcated expressions for the Liouville structure constants 



give 



n 



(x,.7(!f)6'-") 



(Q/6-4) 



T^(6)T^(g/2)T^(g/2 + 6). 



(59) 



The special function Tf,(a;) is the standard element of the LFT (see [19] for the definition 
and properties). Explicitly, 



( 7r//7(— )6 



(Q/b-4) 



exp 



.10 



dt 
T 



3(1 - b^fe-' _ 6sinh'((l-&^)t/(46)) 
2fe2 sinh(t/(26))sinh(6t/2) 



(60) 



where we use "shift" relations (see [19]) for the last upsilon function to extract the additional 
factor (7(^))^. This allows improving the accuracy of the comparison with the analytic 
result, which contains the same factor with a singularity at 6 = 1. The P dependent parts 
are 



ro{P) 



p2T6(6±iP)T6(Q/2±iP) 



Tl{h - g/4 ± iP/2)T2(g/4 ± iP/2)Tl{b + Q/A ± iP/2)Tl{3Q/A ± iP/2) ' 

AP^T,{b ± iP)T,{Q/2 ± iP) 

T2(3V2 - Q/4 ± iP/2)Tl{b/2 + g/4 ± iP/2)T2(l/2/6 + g/4 ± iP/2) ' 



(61) 
(62) 



where we again use the shift relations to move the arguments of all upsilon functions inside 
the strip [0, Q] where the standard integral representation is applicable and we use the 
notation Tft(a; ± y) = Tft(x + y)Tb{x - y) 



ro{P) = P^exp 



dt 



(5 - -Ib'^ + 5//^)2e-* 
262 



(63) 



2cos(f )(cosh((i^) + 2cosh((i^) + cosh(i2^)) - cos(PO(cosh((i^) + 1) - 6 



ri(P) = 4p2exp 



'^dt 
^ 



sinh(t/(26))sinh(6i/2) 
5(1 - b^e-' 
262 + 



(64) 



2cos(f )(3cosh(ii^) +cosh(5^i^)) - cos(Pi)(cosh(^i^) + 1) - 6 

smh{t/{2b))smh{bt/2) 

For the second integral, we analogously find the explicit expressions 

n^L^i^{^-^)b'-''Y^' '\\b'm^-^)^\b'+^-^)b<'-'''^x 



exp 



'^dt 
T 



3(l-62)2e 
262" 



2 ' ' ' 2 
6sinh^((l -62)t/(46)) 

sinh(t/(26))sinh(6t/2) 



(65) 
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and 



ro(P) = P^exp 



dt 



(5 - 186^ + 376^)e-* 
262 



+ (66) 



2cos(f )(cosh(^i^) +2cosh(^i±P^) + cosh(^5^)) -cos(Pt)(cosh(^i^) + 1) - 6 



ri(P) = 4p2exp| 



sinh(t/(26))sinh(6t/2) 
citr (5-2662 + 376^)e-* 

7 [ 2P + ^'^^ 

2 cos(f )(cosh( (g^) + 2 cosh( (^^) + cosh(il^)) - cos(Pt)(cosh((l^) + 1) - 6 

sinh(t/(26))smh(6V2) 

The conformal blocks can be evaluated effectively using a numerical algorithm based on 
the recurrence relations developed in [14, 24-28]. We do not use the elliptic recursion to con- 
struct the necessary conformal blocks here. It turns out that to attain a convincing accuracy 
of the results, we need to know a very few first terms in the g-cxpansion of the conformal 
blocks. This information can be obtained directly using the definition of the conformal blocks 
in terms of the chain vectors (see Appendix A) and the subsequent elliptic transformation. 
Nevertheless, the elliptic representation and especially the form of the prefactor (i.e., the A 
asymptotic of the conformal blocks; see [24]) is very useful. It allows verifying the exphcit 
expressions for the correlation functions by checking the crossing symmetry requirement, and 
in particular, fixing up all the signs, which is difficult to do starting from general principles. 



6 Numerics 

With (55) and (56), calculating reduces to numerically integrating several integrals of the 
general form 

[ \z{l - z)9l{q)Mz)\'d\ (68) 

where J^p{z) is some Liouville conformal block like in (55) or some more complicated com- 
posite expression like in (56). The integrand can be developed as a double power series in q 
and q in accordance with the general expansion 

oo 

z{l-z)et{q):Fp{z) = (16g)"^6,(P)g^ (69) 



r=0 



where a and the coefficients hr are defined by the concrete choice of the function J-p{z). In 
each term, we can integrate in T2 = Im r explicitly with the result in terms of the function 

$(Ar,n= / dV|16g|'^g"g' = , ^/^^'^ , /"'^' cos(7r(r-nx)e-'^^^(2A+r+o^^_ (70) 
Jf 7^(2/1 + r + /) J.^m 
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For (55), we have the sum of six integrals of form (68), and we obtain the series 



I^{b) = 7^7^ ^ (a^I''^ + A^l'"^ + 4''^) + + Af'^ + A^A , (71) 

T—n V / 



L=0 

where in the last sum 

L 



4'^^ = / r,iP)dPj2t'\P)bii{P)^ (^72 + QVS - 1/2, k,L-k), 

k=0 

/•OO ^ 

4'"^^ / n{P)dPj2t''\P)bii{P)^{PV2 + Qy8,KL-k), 

•^0 fe=0 

/"OO 

4''^^= / r2(P)dP5]6i^'^)(P)6K(P)$(PV2 + gV8-l/2,fc,L-^), 

/>oo 

4'°)= / ri(P)dP5^6f'')(P)6i^:ri(P)*(PV2 + gV8,A;,L-A;), (72) 

•^0 fe=0 

poo ^ 

4'^)= / ri(P)dPj]6f'^)(P)6K(P)$(P72 + g78-l,fc,L-fc), 

•^0 fc=0 

/■OO 

4'°^- / r2(P)ciP5^6f''')(P)6i^ri(P)$(P72 + g78-l/2,fc,L-A;). 
•^0 fe=0 

Each term in (72) is suppressed by a factor maxp \qf^, and the series in L converges very 
rapidly in practice. We found that it suffices to sum up to L = 4 to reach the three- to 
four-digit precision (see Table 1). In the case (1,3), we numerically integrate in the same 
way although we now have 18 integrals of form (68). In this case summing up to L = 4 
we were able to reach the two- to three-digit precision (see Table 2). In Figs. 1 and 2, the 
results of numerically evaluating integrals (29) and (31) are shown as circles while the lines 
correspond to the exact result. 
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b 


num. 


E(i'i)(6) exact 


0.999 


0.9959 


0.9960 


0.95 


0.8049 


0.8050 


0.85 


0.4450 


0.4450 


0.80 


0.2799 


0.2800 


1/V2 


0.0001 





0.65 


0.1555 


0.1550 


0.(30 


0.2877 


0.2800 



Table 1: Numerical data for E*^^'^^(6) at = 1. 



b 


num. 


exact 


0.71 


0.0309 


0.0246 


0.69 


-0.1377 


-0.1434 


0.67 


-0.3030 


-0.3066 


0.65 


-0.4623 


-0.4650 


0.63 


-0.6159 


-0.6186 


0.61 


-0.7644 


-0.7674 


0.59 


-0.9096 


-0.9114 


0.57 


-0.9699 


-0.9747 


0.55 


-0.9060 


-0.9075 


0.53 


-0.8409 


-0.8427 


0.51 


-0.7791 


-0.7803 


0.49 


-0.7197 


-0.7203 


0.47 


-0.6621 


-0.6627 


0.45 


-0.6069 


-0.6075 


0.43 


-0.3342 


-0.3282 


0.41 


-0.0384 


-0.0258 


0.39 


0.1827 


0.2622 



Table 2: Numerical data for E^^'^^b) at // = 1. 



15 




0.4 0.5 0.6 0.7 0.8 0.9 

Figure 1: Direct numerical evaluation of reduced integral (29) (circles) versus the exact 
formula (continuous line). 




0.2 0.3 0.4 0.5 

Figure 2: Direct numerical evaluation of reduced integral (31) (circles) versus the exact 
formula (continuous line). 
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Appendix A. Conformal Blocks and the Chain Vectors 

For definiteness, we use the notation for the Liouvillc sector, though the resuhs concerning 
the conformal block are universal (i.e., are independent of the sector). Schematically, the 
contribution of the given conformal family in the four basic OPE can be written as 

N 

N 

V,{z)W,{0) = ^A-A,-A.-i/2 J2 z^'lN),,, (73) 

AT 

N 

where the so-called chain vectors |A^), | A^), |A^), |iV) (with positive integer or half-integer 
N) are the A^th-level descendent contribution of the intermediate state with the confor- 
mal dimension A appearing in the given operator product expansion. The chain vectors 
are completely determined by the superconformal symmetry. Namely, the superconformal 
constraints lead to the recurrence relations 

'g,\N)u = \N^k),,, 

Gk\N),^ = [A + 2kA^-A2 + N-k]\N-k)^2 
for A; > 0. And 

Gk\N)i2 = 1^ - k)i2 + '2A25k,i/2\N - k)u, ^75) 
Gfc|iV),2 = [A + 2kA, - (A2 + 1/2) + N- k]\N^k),^ - 2A24,i/2|iV^),2 

for k > 0. The normalization of the chain vectors chosen in this text is determined by the 
requirements 

|0) = 1, |0j = l, |0j = -l, |0j = (A-Ai-A2). (76) 

Relations (74) and (75) are equivalent to the linear problem for the coefficients determin- 
ing the chain vectors in terms of the Virasoro basis vectors of the same level. These systems 
can be solved numerically up to a rather high level. 
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The necessary s-channel superconformal blocks are defined via the expansions 



T 

•' e,o 

T 

•' e,o 

•F 3,0 
•F 3,0 





0.3 


A 


02 






a\ 


0,3 


A 


a2 


a4 




di 


03 


A 




04 




di 


as 


A 


02 


04 





z \ ~ z 



A-A1-A2 



Ne1,'Z/2 



z \ — z 



z \ = z 



A-A,-A.-l/2 ^ ^^,,(7V|Ar)3,, 

iVez,z/2 
Arez,z/2 



z ] = z 



A-A1-A2-1/2 



5] ;2^12(iV|iV)34. 



wez,z/2 



In the main text, we use the brief notation 



Foo{A,z) = F, 



Fio{A,z)=F, 



F2o{A,z) = F, 



F^o{A,z)=F, 



f a 


a 




A 


V « 


a 




/ d 


d 




A 


V « 


a 




/ a 


a 


A 


V « 


a 




/ d 


a 


A 


V « 


d 





z], Fo^{A,z)=Fo 



f a 


a 




A 


V « 


a 





z], Fi,{A,z)=Fo 



f d 


d 




A 


V « 


a 





z], F2i{A,z)=To 



z], F:,,{A,z)=Fo 



a a 
d a 



a a 
a d 



(77) 
(78) 

(79) 
(80) 



(81) 



(82) 



(83) 



(84) 
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